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; Abstract 

i-G ' We show that an ensemble of ^S molecules in the rotationally ground state trapped on an optical 

lattice exhibits collective spin excitations that can be controlled by applying superimposed electric 
and magnetic fields. In particular, we show that the lowest energy excitation of the molecular 
ensemble at certain combinations of electric and magnetic fields leads to the formation of a magnetic 
Frenkel exciton. The exciton bandwidth can be tuned by varying the electric or magnetic fields. 
We show that the exciton states can be localized by creating vacancies in the optical lattice. The 
localization patterns of the magnetic exciton states are sensitive to the number and distribution of 
vacancies, which can be exploited for engineering many-body entangled spin states. We consider 
the dynamics of magnetic exciton wavepackets and show that the spin excitation transfer between 
molecules in an optical lattice can be accelerated or slowed down by tuning an external magnetic 
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INTRODUCTION 

A major thrust of current experimental research is to create ultracold ensembles of polar 
molecules trapped on an optical lattice [1, 2]. Optical lattices of ultracold molecules are pre- 
dicted to be ideally suited for quantum simulation of complex many-body quantum systems 
[3-5] and the development of new schemes for quantum information storage and processing 
[6, 7]. These applications exploit (i) the possibility to produce molecular ensembles in a 
Mott-insulator state, i.e. an ordered array with one molecule per lattice site [8, 9]; (ii) the 
presence of long-range dipole - dipole interactions that can be used to couple molecules in 
different lattice sites [10]; (iii) the rotational structure of molecules; (iv) the possibility to 
address spectroscopically molecules in specific sites of an optical lattice [6]. For example, 
DeMille proposed to use rotational states of polar molecules confined in a one-dimensional 
optical lattice as qubits of a quantum computer entangled by the dipole - dipole interactions 
[6]. Recently, Micheli and coworkers showed that ultracold molecules in the ^S electronic 
state trapped in a two-dimensional optical lattice can be used for engineering lattice-spin 
models that give rise to topologically ordered states [3]. This stimulated experimental work 
on the creation of ultracold molecules in electronic states with unpaired electrons [1] . 

Ultracold ^S molecules can be produced by photoassociation of ultracold alkali metal 
atoms with ultracold alkaline earth [11] or closed shell lanthanide atoms [12, 13], buffer gas 
loading [14, 15] or direct laser cooling [16], as was recently demonstrated for the molecule 
SrF [17]. The presence of the unpaired electron in a ^S molecule allows for new applications 
of ultracold molecules exploiting weak couplings of the electron spin with the rotational 
angular momentum of the molecule [1-3]. In particular, the electron spin of ^S molecules 
can be used for encoding quantum information as in other spin-1/2 particles [18] and for 
quantum simulation of many-body spin dynamics [3] . In order to realize these applications, 
it is necessary to develop techniques for entangling the spin degrees of freedom and con- 
trolled preparation of many-body spin-dependent states of ultracold molecules on an optical 
lattice. Micheli and coworkers showed that the electric dipole - dipole interaction between 
^S molecules on an optical lattice leads to spin-dependent binary interactions whose pa- 
rameters can be tuned by a combination of dc electric and microwave fields [3]. Here, we 
extend the work of Micheli and coworkers to explore the possibility of tuning the dynamics 
of collective spin excitations in an ensemble of ^S molecules on an optical lattice by external 



electric and magnetic fields. 

We consider SrF(^S) molecules in the ro-vibrationally ground state confined in a Mott- 
insulator state in an optical lattice. In the presence of a weak magnetic field, the lowest 
energy excitation of the molecular crystal corresponds to the spin-down to spin-up transition 
in an isolated molecule. We show that for certain combinations of superimposed electric and 
magnetic fields, the lowest energy excitation of the molecular crystal leads to the formation 
of a magnetic Frenkel exciton. The magnetic exciton is a many-body entangled state of the 
molecules. If some of the molecules are removed from the optical lattice to produce vacancies, 
which can be achieved by a recently demonstrated technique [19], the exciton undergoes 
coherent localization due to scattering by the impurities. Using a series of examples for 
SrF molecules in an optical lattice, we show that the localization of the magnetic exciton 
can be controlled by varying the concentration and distribution of vacancies as well as the 
external magnetic or electric fields. The system proposed here can be used for quantum 
simulation of spin excitation transfer in many-body crystals without phonons. We show 
that the spin excitation transfer between molecules can be accelerated or slowed down by 
tuning an external magnetic field. 

THEORY 

Energy levels of a ^S molecule 

The Hamiltonian of an isolated ^S molecule in the presence of superimposed electric and 
magnetic fields can be written as [20-22] 

H^ = i^ro-vib + IsrS ■N-E-d + 2iibB ■ S, (1) 

where the first term determines the ro-vibrational structure of the molecule, 7sr is the 
constant of the spin-rotation interaction between the rotational angular momentum N and 
the spin angular momentum S of the molecule, E and B are the vectors of the electric 
and magnetic fields, d is the dipole moment of the molecule and fiB is the Bohr magneton. 
We assume that both E and B are directed along the quantization axis z. We consider 
molecules in the vibrationally ground state and use the rigid-rotor approximation for our 
calculations. It is convenient to use the basis of direct products [21] of the rotational |A^Mjv) 



and spin \SMs) wave functions to evaluate the eigenvectors and eigenvalues of Hamiltonian 
(1). Here, Mn and Ms denote the projections of N and S, respectively, on the z axis. The 
diagonalization of Hamiltonian (1) yields the energy levels of the molecule in superimposed 
electric and magnetic fields, shown in Figure 1 for the particular example of SrF. The 
eigenvectors of Hamiltonian (1) are linear combinations 

^f=J2 J2'^nm^Ms\NMn)\SMs) (2) 

NMpj Ms 

The index / denotes the /-th excited state of the molecule. We neglect the hyperfine 
structure of the molecule. This is a good approximation for the magnetic fields considered 
in the present work. 

At zero electric field, the magnitude of the rotational angular momentum N is conserved. 
When N = 0, the spin-rotation interaction vanishes and two lowest energy levels of the 
molecule in a weak magnetic field correspond to the projections Ms = —1/2 and Ms = +1/2 
of the spin angular momentum S on the magnetic field axis. The state with Ms = —1/2 
(state a in Figure 1) is the absolute ground state of the molecule at all magnitudes of the 
magnetic field. The state with Ms = +1/2 (state /3 in Figure 1) becomes degenerate with a 
high-field-seeking Zeeman state (state 7 in Figure 1) of the A^ = 1 manifold at some value 
of a magnetic field. For SrF, this degeneracy occurs at the magnetic field 535.3 mT. States 
P and 7 have different parity. The E ■ d interaction is the only term in Hamiltonian (1) 
that couples states of different parity. Therefore, the crossing between states /3 and 7 is 
real in the absence of an electric field and becomes avoided in the presence of an electric 
field. The properties of ^S molecules, such as alignment and orientation, are very sensitive 
to electric and magnetic fields near these avoided crossings [22-24]. In the present work, 
we show that the avoided crossings depicted in Figure 1 can be exploited for inducing and 
controlling collective spin excitations of ^S molecules on an optical lattice. 

Magnetic excitons induced by electric fields 

We consider an optical lattice of SrF molecules prepared in the absolute (vibrational, 
rotational and Zeeman) ground state with one molecule per lattice site. We assume that 
the lattice sites are separated by 400 nm and the tunneling of molecules between lattice 
sites is suppressed, i.e. the molecules are in a Mott insulator state. This can be achieved 
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FIG. 1: Energy levels of the SrF(^S) molecule as functions of a magnetic field in the presence 
of an electric field of 1 kV/cm. The rotational constant of SrF is 0.251 cm~^, the spin-rotation 
interaction constant 7sr is 2.49 x 10^^ cm^^ and the dipole moment is 3.47 Debye. States /? and 
7 undergo an avoided crossing at the magnetic field value Bq. The value of Bq varies with the 
electric field. 

by applying laser fields of high intensity [9] . With the current technology, it is possible to 
create optical lattices that provide harmonic confinement for ultracold atoms and molecules 
with the vibrational frequencies of the translational motion up to 100 - 150 kHz [25]. We 
assume that the molecules populate the ground state of the lattice potential and neglect 
the center-of-mass motion of the molecules that broadens the spectral lines of the molecular 
crystal by about 5 % [26]. 

The Hamiltonian describing the optical lattice with identical ^S molecules in the presence 
of superimposed electric and magnetic fields is 



A^ 1 N ^ N ^ 

H =^H^^n + -"^"^Vddirn-rm), (3) 

ra=l ra=l rriy^n 

where r„ is the position of the n-th lattice site, Has,n is Hamiltonian (1) for molecule in site 
n, Vdd is the electric dipole-dipole interaction between molecules in different lattice sites, 
and A^moi is the total number of molecules. Hamiltonian (3) can be rewritten in terms of 

\f and Bnj 
follows [27]: 



the operators Bl r and Bn f that create and annihilate a molecular excitation / in site n as 



^ = E^/ + EEEE Ji%Bi^B^j, (4) 



where 



Hf = efY^ BijKf + E E JliBijBmj, (5) 

6/ is the energy of the /-th excited state of the molecule and 

Jn,m = {<t)f{rn)Mrm)\Vdd{rn " n«) |0o(nO0g(^m)) (6) 

The eigenstates of Hamiltonian (5) can be written as 



^/N—^ 



'mol 



where the wave functions 



describe the ensemble of A^moi — 1 molecules in the ground state 0o and one molecule located 
at Vi in the /-th excited state. Eq. (7) gives the wave function for a Frenkel exciton corre- 
sponding to an isolated /-th excitation of the molecule [27]. In general, different excitonic 
states are coupled by the terms Jl'^ with f ^ g in Eq. (4) and the eigenstate of Hamiltonian 
(3) is a linear combination of the excitonic states (7) associated with different excitations /. 
Hamiltonian (5) can be diagonalized by the unitary transformation [27] 

Bl r = -^L= y Bl re-""-"-- (9) 



to yield 

^/ = Ei^/ + ^/(^)]^IA/' (10) 

k 

where k is the exciton wavevector and Lf{k) = ^^ J^'qC*'""^" determines the exciton disper- 
sion in the absence of inter-exciton couphngs. If all off-diagonal couplings J^'f^ with f ^ g 
are much smaller that the energy separation of the molecular state / from other molecu- 
lar states, the coupling terms J^'^ in Eq. (4) can be neglected and Eq. (7) provides an 
accurate description for the wave function of exciton /. This corresponds to the two-level 
approximation [27, 28]. 

In the present work, we consider the lowest energy excitation of the molecular crystal 
corresponding to the Zeeman transition from state a to state /3 depicted in Figure 1. In the 
absence of electric fields and at low magnetic fields B ^ Bq^ state a is |A^ = 0)|M5' = —1/2) 
and state (5 is \N = 0)|M5' = -1-1/2). Because the electric dipole - dipole interaction is 
spin-independent, integral (6) for 0o = 1^ = 0, M^ = —1/2) and ipf = \N = 0,M5 = 
-hl/2) vanishes. The |A^ = 0, M^ = -1/2) -^ \N = Q,Ms = +1/2) excitation therefore 
produces a single isolated molecule in the excited state. Near the avoided crossing shown in 
Figure 1, state /3 is a linear combination of |A^ = 0, M^ = —1/2) and |A^ = 1,M5 = 1/2) 
and the coupling constants J/'^ must be non-zero. The matrix elements J^'4 determine 
the exciton bandwidth defined as the difference Aj = |Lj(7r/a) — Lj(0)|, where a is the 
intermolecular separation. Figure 2 shows the dominant nearest-neighbour coupling constant 
J = Jnn+i ^^^ ^he corresponding exciton bandwidth as functions of the magnetic field for 
three values of the electric field for a one-dimensional array of SrF molecules on an optical 
lattice. As the magnetic field is increased through the avoided crossing, state /3 changes 
from predominantly |A^ = 0, Ms = 1/2) to predominantly |A^ = 1, Ms = —1/2). As a result, 
the coupling constants J^'4 giving rise to the lowest energy exciton as well as the exciton 
bandwidth increase from zero to a finite value. The variation of the coupling constants and 
the bandwidth with the magnetic field is slower for higher electric field. 

The results presented in Figure 1 indicate that the first excited state of the SrF molecule 
is separated from other excited states by more than several MHz. Figure 2 demonstrates 
that the magnitudes of the coupling constants Jl'^ ^-iid the bandwidth of the lowest energy 
exciton are < 20 kHz. The off-diagonal matrix elements J^'^ have similar magnitude as 



Jnm- The results of Figures 1 and 2 thus indicate that the lowest energy exciton is well 
separated from the other molecular transitions so it should be well described by a single 
exciton wave function of the form (7). To prove this we calculated the energy band of the 
lowest energy exciton by the direct diagonalization of Hamiltonian (4) with four coupled 
excitonic states and by the diagonalization of Hamiltonian (5), which corresponds to the 
two-level approximation [27]. Figure 2 demonstrates that the two- level approximation for 
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FIG. 2: Upper panel: Nearest-neighbour coupling constants J = J^n+i (upper panel) giving rise 
to the lowest energy exciton in a crystal of ^S molecules calculated for three values of the electric 
field as functions of the magnetic field. Lower panel: Exciton bandwidths A calculated for three 
values of the electric field as functions of the magnetic field. The calculations are for 2000 SrF(^E) 
molecules arranged in a one-dimensional array with the lattice spacing 400 nm. Lines in the lower 
panel represent the results obtained by the diagonalization of Hamiltonian (4) with four coupled 
excitons. Symbols show the results of the calculation using the two-level approximation for the 
lowest energy excitation. 

Figure 2 suggests an interesting possibility of creating an entangled many-body state 
(7) of non-interacting molecular spins. This can be achieved by the following procedure: (i) 
apply an electric field to couple the molecular states of different parity; (ii) tune the magnetic 
field adiabatically to a value B > Bq (see Figure 1), where the a — )■ /3 exciton bandwidth is 



8 



large; (iii) generate the lowest energy excitation; (iv) detune the magnetic field to a value 
B <^ Bq and turn off the electric field. As mentioned in the previous section, when the 
electric field is absent and the magnetic field is -B ^ -Bqi state a is a pure Ms = —1/2 state 
and state /3 is a pure Ms = +1/2 state. If carried out faster than the inherent time scale of 
the exciton dynamics (see next section), step (iv) must therefore project the excitonic wave 
function on the many-body state 

^ = 7^E^^^f (11) 

with 

^f = \Ms = 1/2),^ n 1^^ = -l/2)-r (12) 

The variation of the magnetic field in step (iv) must be generally faster than h/J to preserve 
the exciton state, but slow enough to preclude the non-adiabatic transitions to molecular 
state 7. This can be achieved because the splitting of the molecular states /3 and 7 at the 
avoided crossing is much greater than J. We have confirmed that the magnetic field can be 
detuned to a value B <^ Bq without changing the magnitudes of the coefficients Cj by time- 
dependent calculations as described in the next section. Wave function (11) is a many-body 
analogue of the two-body EPR state with spin-1/2 particles a\ t)| i) + b\ |)| t) [18, 29]. In 
the following section, we show that the expansion coefficients Cj in Eq. (11) can be modified 
by creating vacancies in the molecular crystal. 

Localization of magnetic excitons in the presence of vacancies 

In an ideal, infinitely large crystal, the exciton wave function is completely delocalized 
and the probability to find any molecule in the excited state is the same. If the molecular 
crystal contains impurities in the form of other molecules or vacancies, the exciton wave func- 
tion is modified and the excitons may undergo coherent localization [27]. The localization 
patterns depend on the dimensionality of the crystal, the concentration and the distribution 
of impurities as well as the exciton - impurity interaction strength. Wurtz and coworkers 
have recently presented the results of an experiment showing that atoms in specific sites of 
an optical lattice can be selectively evaporated by focusing an electron beam onto a par- 
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ticular lattice site [19]. This technique is general and can be applied to molecules as well 
as atoms. This method can be used to create mesoscopic ensembles of ultracold particles 
with arbitrary spatial patterns. In the following two sections, we show that the interactions 
of magnetic excitons (11) with empty lattice sites are very strong and that the entangled 
spin states given by Eq. (11) can be effectively modified by evaporating molecules from the 
optical lattice to create vacancies. In this section, we consider a one-dimensional array of 
2000 SrF molecules on an optical lattice with lattice spacing a = 400 nm. Since the exciton 
coupling constant J is negative (c/. Figure 2), the exciton has a positive effective mass, i.e. 
the energy of the exciton increases with the wavevector [27]. An excitation of molecules 
by laser field produces excitons with small wavevector [27] so we focus on a state near the 
bottom of the energy spectrum. 

To find the eigenstates of a molecular crystal with vacancies, we use the two-level approx- 
imation and modify Hamiltonian (5) and the basis wave functions (12) to omit the terms 
corresponding to the positions of the empty lattice sites. The two-level approximation was 
shown to be very accurate in the previous section. The dimension of the Hamiltonian matrix 
is (A^moi — A^vac) X (A^moi — A^vac), wherc iVvac is the number of vacancies in the lattice. We 
calculate the eigenvectors and eigenvalues using two different numerical methods in order 
to verify the accuracy of the computations. Figure 3 shows the probability density of the 
magnetic exciton for the array of SrF molecules with no defects (panel a), with vacancies 
positioned at every 200-th site (panel b), with a Gaussian distribution of 37 randomly placed 
vacancies (panel c) and with a uniform distribution of 37 randomly placed vacancies (panel 
d). The results show that the presence of vacancies modifies the exciton wave function to 
a great extent and leads to various localization patterns of the exciton wave functions. In 
particular, a periodic placement of vacancies (Figure 2b) leads to the formation of wave 
function domains with the molecules in the middle of the crystal being much more likely 
in the Ms = +1/2 state. The Gaussian and uniform distributions of impurities lead to the 
formation of a strongly localized state with only a few molecules sharing the Ms = +1/2 
excitation. The localization patterns can be further modified by varying the concentration 
of impurities, which is demonstrated in Figure 4. 
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FIG. 3: Magnetic exciton wave function (12) for a one-dimensional array of SrF(^S) molecules in 
an optical lattice with 2000 lattice sites separated by a = 400 nm. Panel (a): Free exciton wave 
function. Panel (b): Exciton wave function in the presence of equally spaced vacancies. Panel (c): 
Exciton wave function in the presence of a Gaussian distribution of randomly placed vacancies. 
Panel (d): Exciton wave function in the presence of a uniform distribution of randomly placed 
vacancies. The number of vacancies for Panels (c) and (d) is 37. The positions of the vacancies in 
Panels (b), (c) and (d) are indicated by the symbols. The wave function is presented for state 51 
in the spectrum of ~ 2000 states. The magnetic field is 536.9 mT and the electric field is 1 kV/cm. 

Dynamics of spin excitation transfer 

The system proposed in the present work may allow for the study of the dynamics of spin 
excitation transfer between molecules in real time. A gradient of a magnetic field applied 
to an optical lattice of molecules initially in the Ms = — 1/2 state can be used to excite a 
molecule in a specific lattice site to the Ms = +1/2 state by an rf tt/2 pulse resonant with 
the Ms = —1/2 — )■ Ms = +1/2 transition for the molecule in this particular position [6]. 
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FIG. 4: Magnetic exciton wave function (12) for a one-dimensional array of SrF(^S) molecules in 
an optical lattice with 2000 lattice sites separated by a = 400 nm in the presence of a uniform 
distribution of randomly placed vacancies. Different panels correspond to different concentrations 
of vacancies. Panel (a): 20 vacancies; panel (b): 40 vacancies; panel (c): 60 vacancies; panel (d): 
80 vacancies. The positions of the vacancies are indicated by the symbols. The magnetic field is 
536.9 mT and the electric field is 1 kV/cm. 

This should produce a spin excitation locahzed on a single molecule. The spin excitation 
transfer can then be induced by removing the field gradient, tuning the magnetic field to a 
value B r^ Bq and applying an electric field that gives rise to significant couplings J^'^ in 
Eq. (5). The spin excitation of a single molecule is a localized exciton wave packet. When 
the coupling matrix elements J^'^ are non-zero, the wavepacket must exhibit dynamics of 
spin excitation transfer determined by the strength of the coupling matrix elements </f '4- 

In order to explore the effects of electric and magnetic fields on the dynamics of the exciton 
wavepackets and determine the timescale of the spin excitation transfer, we consider a small 
ensemble of SrF molecules on a one-dimensional optical lattice. In a finite-size ensemble of 
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molecules, the spin excitation should propagate to the edge of the molecular sample and 
return to the original molecule, leading to revivals of the exciton wavepacket. We expand 
the total wave function of the system in terms of direct products (12) with time-dependent 
expansion coefficients, 

^ = j2Fnmi (13) 

n 

The substitution of this expansion in the time- dependent Schrodinger equation with Hamil- 
tonian (5) yields a system of first-order coupled differential equations for the expansion 
coefficients Fn{t). We integrate these equations to obtain the time evolution of the quan- 
tities |-F„(t)p, which give the probability for the spin excitation to be found at time t in 
lattice site n. Figure 5 shows the autocorrelation function |A(t)P = |(^(i = 0)|^(t))p, 
which describes the probability for the spin excitation to be on the initially excited molecule 
(molecule four in this example) in a system of seven molecules calculated as a function of 
time for four different magnitudes of the magnetic field near the avoided crossing shown in 
Figure 1. The electric and magnetic fields are kept constant in this calculation. Figures 5 
and 6 demonstrate that the spin excitation is transferred to other molecules and then comes 
back to the original site at regular time intervals. The frequency of the revivals is very 
sensitive to the magnetic and electric fields and varies from ~ 1 to ~ 5 ms. 

Vacancies introduce potential barriers which infiuence the dynamics of excitons. The 
strength of the barriers can be modified by producing several vacancies in adjacent lattice 
sites. Figure 7 demonstrates the interaction of the spin excitation with the potential barriers 
produced by two (panel a) and three (panel b) adjacent vacancies. This calculation shows 
that the excitation can tunnel through the vacancies and that the rate for tunneling through 
three empty lattice sites is almost one order of magnitude smaller than that for tunneling 
through two empty sites. Figure 7 demonstrates that the spin excitation in a mesoscopic 
ensemble of ^S molecules in an optical lattice separated by three empty lattice sites can be 
considered isolated on the time scale of < 5 ms. 

SUMMARY 

We have shown that the unique energy-level structure of ^S molecules can be exploited for 
controlled preparation of many-body entangled states (11) of non-interacting molecular spins 
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FIG. 5: Time-dependent probability of the spin excitation to be localized on molecule four in a 
one-dimensional array of seven molecules in an optical lattice with lattice spacing a = 400 nm for 
different magnetic fields near the avoided crossing shown in Figure 1. The electric field magnitude 
is 1 kV/cm. 

on an optical lattice. The proposed procedure is based on inducing the interactions giving 
rise to the formation of magnetic excitons and projecting the exciton states on the basis 
of uncoupled spin states by varying external electric and magnetic fields. Our calculations 
show that the exciton states are formed at particular combinations of electric (Eq) and 
magnetic {Bq) fields near the avoided crossing between states (3 and 7 depicted in Figure 
1. The molecular spin states become uncoupled at magnetic fields B <^ Bq. To confirm 
that the detuning of the magnetic field does not affect the exciton states, we performed 
a time-dependent calculation with seven molecules prepared at t = in the eigenstate of 
Hamiltonian (5) at the magnetic and electric fields near Bq and Eq. The calculation showed 
that the exciton wave function remained the same upon detuning of the magnetic field to 
a value B <^ Bq. This result is largely independent of the speed of the magnetic field 
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FIG. 6: Time-dependent probability of the spin excitation to be localized on molecule four in a 
one-dimensional array of seven molecules in an optical lattice with lattice spacing a = 400 nm for 
different electric fields near the avoided crossing shown in Figure 1. The magnetic field magnitude 
is 536.9 mT. 

variation because the exciton state, once formed, remains the eigenstate of Hamiltonian 
(5) at all values of the magnetic field. It may also be interesting to generate many-body 
entangled states by first creating a spin excitation on an isolated molecule and then tuning 
the magnetic and electric fields close to Bq and Eq. This should produce a dynamically 
changing exciton wavepacket. In order to elucidate the feasibility of projecting the exciton 
wavepackets on the basis of uncoupled spin states by detuning the magnetic and electric fields 
from Bq and Eq, we calculated the time-evolution of the magnetic exciton wavepackets at 
electric and magnetic fields near the values Eq and Bq. Our calculations indicate that the 
time scale for the magnetic excitation transfer at magnetic and electric fields near Eq and 
-Bo is ~ 1 ms. Our results show that it is sufficient to change the electric field by 0.1 kV/cm 
and the magnetic field by 5 Gauss to suppress significantly the time evolution of the exciton 
wavepackets. This variation of the dc magnetic and electric fields can be achieved on the 
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FIG. 7: Time-dependent probability of the spin excitation to be localized on molecule five (solid 
blue line) and molecule 6 (solid red line) in a one-dimensional array of ten molecules in an optical 
lattice with lattice spacing a = 400 nm. Panel (a): molecules five and six are separated by two 
empty lattice sites. Panel (b): molecules five and six are separated by three empty sites. The 
magnetic field magnitude is 536.9 mT and the electric field is 1 kV/cm. At t = 0, the spin 
excitation is localized on molecule four. 

time scale of ~ 1 /is [30, 31]. 

We use SrF molecules for our illustrative calculations in this paper. SrF has a relatively 
large dipole moment, a relatively small rotational constant and a very weak spin-rotation 
interaction by comparison with other ^S molecules [32]. The effects studied here should be 
more pronounced, making the experiments easier, in an ensemble of molecules with a larger 
dipole moment and a larger spin-rotation interaction constant or with an optical lattice 
with smaller lattice-site separations. The experimental work on the creation of ultracold 
^E molecules that should exhibit similar behavior as SrF is currently underway in several 
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laboratories [15-17, 33]. We note that molecules in the ^E electronic state exhibit similar 
avoided crossings in combined electric and magnetic fields [34, 35], which significantly widens 
the range of molecules that can be used for the experiments proposed here. 

The absence of interactions between the spin states in the coherent superpositions (11) 
should make the entanglement robust against decoherence due to vibrational motion of 
molecules in the lattice. We note that this entanglement is created over distances spanning 
thousands of molecules in an optical lattice separated by 400 nm, which amounts to mm size 
scales. We showed that the exciton states can be localized by creating vacancies in the optical 
lattice. Our results demonstrate that the localization patterns of the magnetic exciton 
states are sensitive to the number and distribution of vacancies. This can be exploited 
for engineering entangled states with different patterns. Excitons determine the optical 
properties of crystals and it might also be interesting to explore the interaction of the 
magnetic excitations with electromagnetic radiation, which can be achieved either through 
coupling to a two-photon optical field [36] or to a single photon rf field [37]. The localization 
of excitons by vacancies in the presence of tunable dc electric and magnetic fields can then 
be used for controlled preparation of polaritons. The experimental techniques for producing 
optical lattices of ultracold molecules in the ro-vibrationally ground states [9] and for creating 
vacancies in optical lattices with controlled spatial arrangements [19] have been recently 
demonstrated. 

The system proposed here can be used to study dynamics of spin excitation transfer in 
molecular aggregates. We have shown that localized magnetic excitations in finite-size arrays 
of molecules exhibit revivals in the presence of electric and magnetic fields near Eq and Bq. 
The frequency of the revivals can be controlled by varying the magnitude of the electric 
and magnetic fields. Our calculations show that the magnetic excitation can tunnel through 
multiple vacancies in adjacent lattice sites. The tunneling rate is dramatically suppressed 
when the number of vacancies in adjacent lattice sites is increased. Our results indicate 
that the time scale for the magnetic excitation transfer through three adjacent vacancies is 
> 5 ms. This suggests a possibility of studying spin excitation transfer by measuring the 
exciton wavepacket revivals in small molecular aggregates. The localized spin excitations in 
an optical lattice of ^S molecules can be created and observed by applying a gradient of a 
magnetic field and performing a spatially-resolved spectroscopic measurement [6, 38]. 
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